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We study the gravitational field equations in the presence of a coupling between the derivative of 
a massless scalar field and the Einstein tensor. This configuration is motivated by Galileon gravity 
as it preserves shift invariance in the scalar sector. We analytically obtain solutions with static 
and spherically symmetric geometry, which also include black holes with a single regular horizon. 
We examine the thermodynamical properties of these solutions, and we reveal the non-perturbative 
nature of the coupling constant. We also find a phase transition, similar to the one described by 
Hawking and Page, which occurs at a critical temperature determined by both the black hole mass 
and by the strength of the coupling. 
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I. INTRODUCTION 

In recent years, many extensions of general relativity 
have been considered in the attempt to explain dark en- 
ergy. In particular, a lot of work has been devoted to the 
most general tensor-scalar action that generates equa- 
tions of motion with second-order derivatives, discovered 
many years ago by Horndeski In fact, this action 
shows interesting self-tuning cosmological solutions [2, 
and includes Galileon gravity @ and massive gravity y| . 
Moreover, it can be obtained by compactifying a suit- 
able truncation of higher-dimensional Lovelock gravity 
Q. Among the many interesting features of Horndeski 
gravity there is the coupling between the derivative of 
a scalar field and the Einstein tensor. On cosmological 
backgrounds, this term leads to accelerated expansion 
without the need of any scalar potential, as noted for the 
first time in Q . A similar accelerating effect was also ob- 
served in the case when the coupling involves the Ricci 
tensor alone Q. These features attracted much interest 
in both inflationary pi-flfj| and late-time cosmology fTTj ) . 
but little attention has been paid to local and stationary 
solutions so far. 

In this paper, we wish to partially fill this gap by con- 
sidering static geometries with spherical symmetry in or- 
der to find exact black hole solutions with regular hori- 
zons. Approximate solutions in the presence of an elec- 
tromagnetic field were found in 12]. In a recent paper, 
it has been shown that a no-hair theorem is at work for 
Galileon gravity. In other words, there cannot be asymp- 
totically flat black hole solutions unless the scalar field 
is globally constant Based on this result, Nicolis 

and Hui have proposed an interesting test to constrain 
the parameters of Galileon gravity [L4|]. To evade the 
no-hair theorem we need to relax at least one hypothesis 
and the most natural one is asymptotic flatness, an idea 
supported by the existence of asymptotically de Sitter 
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solutions. 

This strategy turns out to be successful and we are able 
to find analytic solutions with spherical symmetry and a 
regular horizon to the equations of motion. The paper is 
organized as follows. In Sec. 2 we obtain the equations of 
motion and we show that spherically symmetric solutions 
with a regular horizon exist in an analytic form. In Sec. 

3 we focus on the the thermodynamical properties of the 
black hole and the rle of the coupling parameter. In Sec. 

4 we briefly discuss other solutions in different areas of 
the parameter space and we conclude in Sec. 5 with some 
remarks. 



II. EQUATIONS OF MOTION AND ANALYTIC 
SOLUTIONS 

To begin with, let us consider the Lagrangian 
m 2 1 / 7 \ 

where m p is the Planck mass, z a real number, G the 
Einstein tensor, p a scalar field, and g^ v is the metric, 
chosen with mostly plus signature. The absence of scalar 
potential allows for the shift symmetry p — > (^+const, 
which is the relevant Galileon symmetry that survives in 
curved space [j| . For this reason we will refer to (p as to 
the Galileon field in the following. In the usual Galileon 
terminology, the three terms in L are representative of 
L2, L4 and L5, see e.g [la ]. The term involving the Ein- 
stein tensor also appears in the context of massive grav- 
ity, where the parameter z is related to the inverse of the 
graviton mass Our goal is to obtain static solutions 
with spherical symmetry, so the metric ansatz is 

ds 2 = ~F{r)dt 2 + G{r)dr 2 + p 2 {r){d6 2 + sm 2 6d(/) 2 ). (2) 

To obtain the equations of motion, we find more con- 
venient to express the Lagrangian in terms of the metric 
components and then vary the action S = J d^x^fgh with 
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respect to the four fields <p, F, G, p. Once the equations 
of motion are known, one can set p — r and find 



F' 



G-l 



mlr 2 G KralG 2 



F : ztpVFG' 

F' _ 2m p G(G - 1) + zi> 2 (G - 3) + m 2 p r 2 Gip 2 



F 



£1 

F 



G 



(3z</> 2 + 2m%G) 
2m p G(G - 1) - 2zV 2 



zr(4> 2 )' 



(3z-ifj 2 + 2m\G) 



(3) 
(4) 
(5) 



where K is an integration constant and if> = ip'. There is 
also a fourth, redundant equation whose form is unimpor- 
tant. We immediately note that ip = implies K = and 
the resulting metric turns out to be the Schwarzschild 
one. If z = (and tp ^ 0), one finds the Just solution, 
which is known to be singular both at the origin and at 
the horizon, see e.g. [T(|. Finally, we note that the first 
equation includes the term y/FG. Therefore, for K ^ 0, 
G and F must have the same sign for all r for the metric 
components to be real-valued. 

When K = and z ^ 0, we can analytically find exact 
solutions to the system ([5]). Their form crucially depends 
on the sign of z and in this letter we mainly focus on the 
case z > as it is the most interesting. We will comment 
on the case with negative z at the end. The solution to 
the system reads 
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where we defined I 2 = 12z/m 2 and M is a constant of 
integration that will play the role of a mass. There is a 
second constant that multiplies F(r) and that can be ab- 
sorbed into a redefinition of the time coordinate. We im- 
mediately note that the function F is very similar to the 
gu component of a Schwarzschild Anti-de Sitter (SAdS) 
black hole with spherical horizon pj], [HI ■ The analysis 
of the curvature invariants reveals that these are all fi- 
nite for r > 0, and, in particular, at r — r^, namely at 
the zero of the function F(r), which is unique if M > 0. 
On the opposite, the Ricci scalar diverges at r = con- 
firming that there is a physical singularity at the origin. 
Thus, the solution above describes a genuine black hole 
with one regular horizon located at r — if AI > 0. In 
contrast, when M = the metric is non-singular for all 
r > and this vacuum solution will be very important 
for the analysis of the thermodynamical properties [24l |. 

From Eq. ([5]) we see that ip vanishes only when z — » oo. 
In fact, in this limit we recover the Schwarzschild solu- 
tion with F = G- 1 = 1 - 2M/(m 2 p r) and = 0. This 
is also due to the fact that the Galileon term is constant 



on shell as (g rr — zG rr /m 2 )i/; 



-rrip/z for all r. All 



these elements indicate that z is a non-perturbative pa- 
rameter when we regard the Lagrangian (|TJ) as a theory 
of modified gravity. Indeed, the deviation from general 
relativity vanishes when z diverges and the Galileon field 
is strongly coupled. Also, the parameter z clearly in- 
terpolates between the flat black hole solution and the 
SAdS one asl/z essentially plays the role of an effective 
negative cosmological constant. 



III. THERMODYNAMICAL PROPERTIES 

The similarity between the geometry of our solution 
and the one of the SAdS black hole suggests standard 
techniques to study the thermodynamical properties [l8| . 
First of all, the inverse temperature ft is determined by 
the periodicity of the Euclidean metric obtained by the 
analytic continuation t — > —ir, that is [l9| 
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For comparison, we recall that the inverse temperature 
for a SAdS black hole with spherical horizon is given 
by p = 4:TTl 2 r h /(l 2 + 3r 2 h ) [H|. With the definition of 
I given above, one sees that, in the large m p rh limit, 
or for small z, the two /3's coincide. On the opposite, 
in the large z limit we recover the inverse temperature 
associated to the Schwarzschild black hole. For fixed z, 
and similarly to the SAdS case, the temperature diverges 
for M — ¥ and reaches the absolute minimum T m i n = 
m p /(27rv / 2z) — 0.1125 m p /y/z. For large mass, it grows 
again and linearly with r%. 

As for the SAdS black holes, the volume of the Eu- 
clidean action, obtained by the analytic continuation 
t — > —ir, is infinite. Therefore, the partition function 
related to the volume V and to the Helmotz free energy 
A by the relations \nZ = —V = —f3A is meaningless. 
However, along the lines of [18(, we can define a finite 
quantity by subtracting the volume of the action with 
mass M = to the one with M ^ according to the 
formula 



(10) 



V = V o [0o L{r)dr - /3 / L{r)dr 



where Vo is the volume of the horizon space and f> r^. 
The constant (3q is identified with the arbitrary inverse 
temperature that corresponds to the M = background. 
It is related to j3 by requiring that the periodicity of g tt 
for the two backgrounds is the same for any large f, i.e. 
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After some calculations, we find that, in the limit r — > oo, 
the normalized volume reads 
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where we have defined the dimensionless variable x = 
mpfh/ \fz- From this result we can compute the energy 
E = dV/d/3, the entropy S = (3E — V and the heat 
capacity C = dE/dT. By utilizing the implicit relation 
between the mass and the horizon radius of the black 
hole 
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It is easy to show that all these expressions tend to the 
values associated to a Schwarzschild black hole when z 
diverges. We also note that the z — > (z — > oo) limit is 
the same as the ru — > oo (r/j — > 0) one, denoting a sort 
of duality between z and (and hence M). For large z, 
the leading term of the entropy is equal to 7rm 2 r 2 namely 
a quarter of the horizon area. For z —> the value is just 
one half. Most importantly, for finite values of z the en- 
tropy does not follow the area rule and this is reminiscent 
of the effects on the entropy induced by high order cor- 
rections to the Einstein-Hilbert action, such as in Gauss- 
Bonnet gravity (23). This fact also reinforces the con- 
nection between Galileon theory and higher-dimensional 
truncated Lovelock theory discussed in Q. The en- 
ergy E and the entropy S are always positive, except 
in the interval defined by 0.979 \fz < m, p rh < v2z and 
0.883 y/z < m p rh < y/2z respectively. 

The heat capacity C becomes negative for large z and 
tends to the Schwarzschild value — 27rm 2 r 2 . Instead C 
tends to the positive value 47rm 2 r 2 for vanishing z. As 
in the SAdS case, it diverges at the minimum temper- 
ature of the black hole T m i n . The sign of the heat ca- 
pacity and the behavior of the temperature are shown 
in Fig. ([]}. We see that for every T > T min there are 
two black hole solutions with the same temperature. For 
T > T2 the small- x black hole has negative heat capac- 
ity, while the large- a; one has positive heat capacity. We 
recall that small x means small black hole mass or large 
z, that is when the solution tends to the Schwarzschild 
metric. Thus, for T > T2 the situation is the same as for 
the SAdS case: the small- a; black hole decays either into 
radiation or into the large- a: dual [llj]. Between T\ and 
T2, however, there is a region where both black holes are 
thermodynamically stable as C > 0, and this marks a 
difference with respect to the SAdS case. Moreover, be- 
tween T m i n and Ti, the situation is reversed, as the small- 
x black is stable while the large- a; is not. The Helmholtz 
free energy A = j3~ l V is positive only in the tiny interval 
T m tn < T < T h , where T h ~ 0.1131 m p /^/z. 



In the SAdS case, for T min < T < T\ the black hole 
with the largest mass is locally stable. However, the evap- 
oration reduces the mass and the temperature until T m i n 
is reached and, as the free energy is positive here, the 
phase transition to pure radiation occurs (l8l | . The other 
black hole solution, with the smallest mass, is unstable in 
any case as the heat capacity is negative. Therefore, low 
temperature SAdS black holes are naturally attracted to- 
wards the phase transition. In our case the situation is 
more involved as, in opposition to the SAdS case, the 
small black hole is stable (C > 0) in the regime where 
the temperature increases when x decreases, eventually 
keeping the system outside the region where A > 0. It 
is therefore possible that the phase transition to pure 
radiation does not occur at this stage. However, as x 
keeps decreasing, the black hole temperature eventually 
crosses the value T2 at which the heat capacity becomes 
negative. Here, the black hole can dissolve into pure ra- 
diation, which collapses into a large x black hole with the 
same temperature just as in the SAdS case. Therefore, at 
the end of the the day, the black hole will always end its 
thermodynamical cycle on the right branch of the tem- 
perature curve, where x and T decrease simultaneously, 
bringing the system towards the usual phase transition 

at Eynj.ri. • 
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FIG. 1: Plot of the temperature in units of m p /y/~z versus 
x — m v rh/ ' \fz. In the grey regions the heat capacity is neg- 
ative. The tiny horizontal strip just above T m i„ represents 
the region where the free energy A is positive. The relevant 
quantities are: T min = 0.1125 m p /y/z, Ti = 0.118 m v l\fz 
T 2 = 0.147 m, 



xi = 0.655, x div = V2,x 2 = 1.904 
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IV. OTHER SOLUTIONS 

To conclude the analysis with z > 0, we discuss the 
properties of the scalar field. It is easy to see that tjj 2 
is negative for all r > r/ l; which implies that ip is imagi- 
nary. The question is whether this implies instabilities in 
our solution. Examples of black holes with scalar fields 
that become complex in some radial range are know in 
low-energy string theory, see e. g. [2l|. Generally speak- 
ing, in these models the scalar field is complex inside the 
event horizon, while in our case it is the opposite, so 
one might worry that the solution is unstable outside the 
event horizon. Our point of view is that tp should not be 
considered as a matter field but rather as an extra de- 
gree of freedom, expressed by the real quantity i/> 2 . This 
is evident also from the equations of motion ([5]), where 
the new degree of freedom appears only as ip 2 . 

For completeness, we briefly comment on the solu- 
tion with negative z. The function F is basically the 
same as in Eq. (j6]) with the arctan term replaced by 
arctanh(TO p r/y / z), which reduces the domain of r in- 
side the interval [0, y^z/m p ]. One can find a solution 
that extends outside this domain by using the formula 
arctanh(a;) = l/21n[(l + x)/(l — x)] and then allowing 
1 — x — > \1 — x\, where x — m p r/yfz as usual. The result 
is that, in the large r limit, the solution tends to a de 
Sitter space expressed in static coordinates, as expected. 
However, the absolute value in F implies that the met- 
ric is non-differentiable at i = 1, which, in principle, 
calls for a non-trivial stress tensor to be added to the 
Lagrangian. We leave to future work a detailed analysis 
of this solution. 

We conclude this paper by discussing the solutions to 
the system ((5} with K ^ 0. With a simple argument 
we can show that there cannot be black hole solutions in 
this case. The equation of motion ([3]) can be written as 
a conservation of a current, d r J r = 0, where 
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m 2 G 2 



-K. 



(17) 



Now, a regular horizon located at r — implies F(rh) = 
0. But then J r must necessarily vanish for all r, which 
means that K must vanish too. Therefore, black hole 
solutions are possible only when K vanishes. However, 



there could be other interesting solutions. We first look 
at Ricci flat solutions, for which F = G _1 = 1 + Q/r 
for some constant Q. Direct substitution into the system 
([5]) shows that the only consistent solution requires again 
K = 0. Finally, by using numerical integration, we find 
a metric that is regular everywhere except at the origin 
Unfortunately, it turns out that some curvature invari- 
ants diverges precisely at r = so we conclude that when 
K ^ all spherically symmetric solutions host a naked 
singularity at the origin. Interestingly, a similar behavior 
was also observed in the context of intersecting branes in 
higher-dimensional gravity j23| . 

V. CONCLUSIONS 

In this paper we have explored a brand-new class of 
black holes that are exact vacuum solutions of scalar- 
tensor gravity with non-derivative coupling, which is a 
typical feature of Galileon gravity. They show a locally 
asymptotic anti-de Sitter geometry and a rich thermody- 
namical structure, with multiple phase transitions that 
depend on the mass and on the coupling parameter. The 
latter turns out to be non-perturbative, in the sense that 
the usual Schwarzschild solution is recovered when it di- 
verges. 

We believe that our results deserve and need future 
developments. For example, the classical stability of the 
solution should be studied through the analysis of the 
(coupled) fluctuations of the metric and of the Galileon 
field. From an astrophysical point of view, it would be in- 
teresting to choose a coupling between matter fields and 
the Galileon to study the Tolman-Oppenheimer-Volkoff 
equation and the stability of stars. The solutions that we 
have found can be also used to test the theory against so- 
lar system experiments, in order to constrain the parame- 
ter z. The asymptotic structure of the solution might be 
interesting in the context of AdS/CFT correspondence 
and in brane cosmology, along the lines of (22j. Finally, 
black hole solutions should be sought also when more 
terms belonging to the Galileon family are present in the 
Lagrangian, such as powers of Dip. We hope to report 
soon on these issues. 
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